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■ Abstract 



Shape invariance is an important ingredient of many exactly solvable quantum me- 
chanics. Several examples of shape invariant "discrete quantum mechanical systems" 
are introduced and discussed in some detail. They arise in the problem of describ- 
ing the equilibrium positions of Ruijsenaars-Schneider type systems, which are "dis- 
crete" counterparts of Calogero and Sutherland systems, the celebrated exactly solv- 
able multi-particle dynamics. Deformed Hermite and Laguerre polynomials are the 
typical examples of the eigenfunctions of the above shape invariant discrete quantum 
mechanical systems. 

o 



1 Introduction 



Many exactly solvable quantum mechanical systems, for example, the harmonic oscillator 
without/with a centrifugal potential, the coulomb problem, the trigonometric and hyper- 
bolic Poschl- Teller potential, the symmetric top etc. are shape invariant!^. Here we will 
show that this concept is also useful in a wider context of "discrete quantum mechanics" , 
in which the momentum operator p = —ihd/dx appears as an exponential (hyperbolic) 
function instead of a polynomial in ordinary quantum mechanics. For demonstration, 
we present several explicit examples of shape invariant discrete quantum mechanics in 
some detail. The corresponding eigenfunctions are orthogonal polynomials belonging to 
the family of Askey-scheme of hypergeometric orthogonal polynomials [21 El- The exam- 
ples, roughly speaking, are related to one and two parameter deformation of the Hermite 
polynomials and two and three parameter deformation of the Laguerre polynomials. They 
could be interpreted as describing deformed harmonic oscillator with (Laguerre) and with- 
out (Hermite) the centrifugal potential. They also arise in many contexts of theoretical 
physics [Uini. 

These examples arise within the context of multi-particle exactly solvable quantum 
mechanical systems, in particular, the Calogero and Sutherland systems [H] and their 
deformation called the Ruijsenaars-Schneider- van Diejen systems jj]. Their integrability 

1 Contribution to a special issue of Journal of Nonlinear Mathematical Physics in honour of Francesco 
Calogero on the occasion of his seventieth birthday. 
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is based on the root systems and the associated Weyl (Coxeter) group symmetry. More 
than twenty years ago, Calogero 8 found out that the equilibrium position of the A-type 
Calogero system was described by the zeros of the Hermite polynomial. The same fact 
was discussed by Stieltjes in a slightly different context more than a century ago 0. The 
equilibria of the B, BC and D type Calogero systems are described by the zeros of the 
Laguerre polynomials. The equilibria of the Sutherland (trigonometric potential) systems 
based on the classical root systems are described by the zeros of the Chebyshev and 
Jacobi polynomials. The interesting dynamics associated with the equilibria of exactly 
solvable multi-particle systems and the polynomials describing the equilibria of Calogero- 
Sutherland (C-S) systems based on the exceptional root systems are presented in [TO] . 
Ruijsenaars-Schneider (R-S) and van Diejen introduced integrable deformation of the C-S 
systems, with several additional parameters The equilibrium positions of the R-S and 
van Diejen systems are discussed by Ragnisco-Sasaki and Odake-Sasaki As expected 
certain multi-parameter deformation of the Hermite and Laguerre polynomials describe 
the equilibrium positions of the rational R-S or the van Diejen systems. 

The single particle dynamics related to these polynomials is the main subject of this pa- 
per. They are certain deformation of the harmonic oscillator without/with the centrifugal 
potential belonging to the "discrete quantum mechanics" . We show that they are exactly 
solvable thanks to the inherited shape invariance. To the best of our knowledge, they are 
the first examples of shape invariance discussed in the framework of discrete quantum me- 
chanics. Of course there are plenty of works applying the factorisation method |12M13llT4"| 
to various difference equations |15| I16j. However, all of these difference equations contain 
finite shifts in the real direction, which result in polynomials of a discrete variable, for 
example the Charlier, the Meixner and the Hahn polynomials |2j. Apparently they look 
very different from those occurring in quantum mechanics. 

This paper is organised as follows. In section two, the simplest example of a shape 
invariant discrete quantum mechanics, related to a deformed harmonic oscillator is intro- 
duced and discussed in some detail. The eigenfunctions are deformed Hermite polynomials, 
or a special case of the Meixner- Pollaczek polynomials J3J. The course of the arguments is 
essentially the same as in ordinary quantum mechanics or the Sturm-Liouville problem, as 
most clearly formulated in the seminal work of Crum ^3] • It is also known as the factori- 
sation method |12j . or the supersymmetric quantum mechanics |14j . Here we follow the 
notion and notation of Crum's paper. In section three we present three explicit examples 
of shape invariant discrete quantum mechanics, in which a two parameter deformation of 
the Hermite polynomials, a two and three parameter deformation of Laguerre polynomials 
are the eigenfunctions. They are a special case of the continuous Hahn polynomial, the 
continuous dual Hahn polynomial and the Wilson polynomial [Hj- Shape invariance is 
demonstrated elementarily. The meaning of the factorisation at the level of the Hamilto- 
nian and at the level of the polynomial solutions are elaborated in detail. Although these 
polynomials are perfectly well understood in their own right, we do believe dynamical 
interpretation in terms of the factorised Hamiltonian, the analogues of the creation and 
the annihilation operators, etc. would shed new light on them. Section four is devoted to 
the clarification of the background of the present research: how these polynomials came 
to our notice through the polynomials describing the equilibrium points of the van Diejen 
systems based on the classical root systems. The final section is for comments, a summary 
and an outlook. 
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2 Deformed Hermite (Meixner-Pollaczek) Polynomials 

The simplest and best known example of exactly solvable and shape invariant quantum 
mechanics is the harmonic oscillator 

H = -\dl + \(x 2 -l)= l - (-id x + ix) {-id x - ix) , d x = A (2.1) 

The Hamiltonian is factorised and the eigenfunctions are the Hermite polynomials, H n (x): 



H(j) n = £ n (j) n , £ n = n, 4> n (x) oc H n (x) e x ' , n = 0, 1, . . . , . (2.2) 

One naively expects that a certain deformation of the Hermite polynomials would 
constitute the eigenfunctions of a simplest shape invariant "discrete quantum mechanics" . 
This is exactly the case and will be explained in some detail in this section. 

Let us start with the following Hamiltonian, 

H = ^^7^) + ly/V^e^y/Vfx)- \{V{x) + V*(x)), (2.3) 
V(x) = X + ix, V*{x) = X-ix; X G M+. (2.4) 

In the discrete quantum mechanics the momentum operator p = —ihd x (with H = 1) 
appears as exponentiated instead of powers in ordinary quantum mechanics. Thus they 
cause a finite shift of the wavefunction in the imaginary direction. For example, 

e ±i9x 0(x) = <f>{x±i). 

Throughout this paper we adopt the following convention of a complex conjugate function: 
for an arbitrary function f{x) = X] n a n £ n > a n £ C we define f*(x) = Yl n a n xn - Here c* 
is the complex conjugation of a number c G C. Note that f*(x) is not the complex 
conjugation of f(x), (f(x))* = f*(x*). This is particularly important when a function is 
shifted in the imaginary direction. 



Factorised Hamiltonian It is easy to see that the above Hamiltonian ()2.3|) is factorised: 
H = A^A, (2.5) 
A = A(x;X) = ^(V^y^^-e^V^O), ( 2 -6) 

def 1 



A* = A{x;\y = ^(yV{x)e-^ - y/V^e^ 8 "). (2.7) 

Here f denotes the ordinary hermitian conjugation with respect to the ordinary L 2 inner 
product: (/,<?) = (f(x))*g(x)dx. Obviously the Hamiltonian ()2.3|) is hermitian (self- 
conjugate) and positive semi-definite. 

The ground state 4>q is annihilated by A: 



Act) Q {x) = 0, 4> {x) oc y / r(A + ix)T(X - ix) =^ H^> [x) =0, £ = 0. (2i 
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The above equation reads 



'V*(x - §) - I) = \JV(x + §) + §)• (2-9) 

The other eigenfunctions of the Hamiltonian (|2.3j) can be obtained in the form 

Hcj) n = £ n (j) n , <p n (x) oc P n (x)(j)o(x) HP n = £ n P n , (2.10) 

in which if is a similarity transformed Hamiltonian in terms of the ground state wave- 
function (j)Q : 

H = f 0q 1 o H o cf)Q 

= \V{x)e- id * + \V*{x)e id * - \{V(x) + V*{x)). (2.11) 

The eigenvalue equation (|2.10|) is now the following difference equation 

(A + ix)P n (x - i) + (A - ix)P n (x + i) = 2(A + £n)Pn(x), (2.12) 

which obviously has a polynomial solution, P n (x) = a n x n + lower degree terms, of definite 
parity, P n (—x) = (—l) n P n (x). By comparing the coefficient of the highest degree x n term, 
we find easily 

£n = n, n = l,2,...,. (2.13) 
The orthogonal polynomials {P n (x),n > 0} with respect to the weight function 

u(x; A) = T(A + ix)T(X - ix) oc o (x) 2 (2.14) 
satisfy the three term recurrence, which reads 

p n +i(x) = xp n (x) - n(n + 2A - l)/4p n _i(x), n > 0, po = I, p~i = 0, (2.15) 
for the monic polynomial, p n (x) = x n + lower degree terms. With proper normalisation 

2™ 7T 

Pn{x) = -r Pn (x) = (x; -), (2.16) 

n\ I 

it is a special case of the Meixner-Pollaczek polynomial. Here we follow the notation of 
Koekoek and Swarttouw pJJ. Since the second argument ^ remains fixed throughout our 
discussion, we will omit it: 

^ A) (x;f) d = f Pi A) (x). 

They are polynomials in both x and A with real and rational coefficients. The Meixner- 
Pollaczek polynomial also appears in a "relativistic oscillator" model 
Corresponding to the factorisation of H (|2.5[) . H is also factorised: 

H = BC, (2.17) 

C = C(x) = §(e-2 & -ei 9x ), (2.18) 

B = B{x; A) = -i(V(x)e'i 9x - V*(x)e^ dx ) = -i ((A + ix)e~^ - (A - ix)e^ . 

(2.19) 
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They factorise the eigenvalue equation 

BCP£\x) = nP,W(x) <= CPW(x) = P^^x), PP^W = ^ A) (*), (2-20) 
giving rise to the forward and backward shift relations. They read explicitly: 

C: pW(x + i)-Pi A )(x-§) = 2iP^ ) (x), (2.21) 
P: i(^X-ix)P^\x + ^-(X + ix)P^\x-l^=nPi x \x). (2.22) 

Let us define a new set of wavefunctions 

ltn = Acj) n , n = l,2,...,. (2.23) 
As a consequence of the factorisation, they form eigenfunctions of a new Hamiltonian 

Pi = AA^ (2.24) 
with the same eigenvalues {£ n }'- 

Pi0i,„ = AA^A<j) n = A£ n <f> n = £ n (j)i,n, n = 1, 2, . . . , . (2.25) 

Shape Invariance It is straightforward to evaluate the reversed order product AA^: 

Pi = \ e~ id * sjv*(x) + \ yjv*(x) e id * ^/VxJx) - ± (Vx (x) + V? (x) ) + £i , 

(2.26) 

Vi(a;)=A+|+tx, ^*(x) = A + ± - ix, £x = 1. (2.27) 

It has the same /orm as P with a constant term £\ = 1 added and the parameter A shifted 
by \. Thus it is again f actor ised: 

H x = AA^ = A\A x + £x, (2.28) 
A! = Ax(x; A) = {e- L ^ Jv*(x) - e& s/VxJx)) = A(x; A + ±), (2.29) 

4 = Ai(x;A)t d = -^=(7^(x7e-i & - y /v 1 *(x)e^ = A\x;\ + ±). (2.30) 

The ground state wavefunction is ^i annihilated by Ai: 



Ai^ 1) i(z) = 0, (x) oc yr(A + i + ix)r(A + i -ix) (2.31) 
=> Hx4>x,i{x) = £x<f>i,i{x), Ex = 1. (2.32) 

Clearly this process can be repeated as many times as the number of discrete levels of 
the original Hamiltonian P (|2.3|) : 

H = H ^H 1 ^H 2 ► #,->.... (2.33) 
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All these Hamiltonians share the same spectra {£ n }, n = 0, 1, . . .. The eigenfunction of 
the s-th Hamiltonian H s are given by 

def 

H s (j) Sin = £ n ^s,n, <^s,n = A s _iA s _ 2 ■ ■ ■ A\A <f) n , < s < n, (2.34) 
in which 

H s = A._ 1 Aj_ 1 + £ s -i = A\A S + £ s , (2.35) 

A s = A s (x; A) = (e^ 9 * - ^/VJx)) = A(x; A + §), (2.36) 

At = A s (x;X^ d ^ i(V^)e-^-^e^) = A(:r;A+§)t, (2.37) 

F s (x) = A + § + V7(jc)=A + f-ix, £ s = s, (2.38) 

A a fl , s = oc ^/r(A+ f + ix)r(A + § -ix) (2.39) 

=>- H s (f> SjS (x) = £ s (f) SjS (x), £ s = s. (2.40) 

By multiplying Al to s +i, n = f A s <fi S)n we obtain 

Aj0 s +i,n = A\A s (j) s ^ = (H s - £ s )4>s,n = (£n ~ £ S )<t>s,n- ( 2 -41) 



Since the ground state (j) njn of the n-th Hamiltonian H n is known explicitly, we can express 
the n-th eigenfunction (j) n = 0o in of the original Hamiltonian H (|2.3|) in terms of A^ and 
§n,n by repeated use of the above formula (|2.41|l : 

M*; A) = -4- • • • A }-} ^ n (x; A) = ±A\a\ • • • A^M^ A+f )• (2.42) 

This is another formula giving the eigenfunction (j) n (x; A) oc Pn ' (x)(/)q(x; A). The situation 
is depicted in Fig.l. The operator A acts to the right and A^ to the left along the 
horizontal (isospectral) line. They should not be confused with the annihilation and 
creation operators, which act along the vertical line of a given Hamiltonian H s going from 
one energy level n to another n db 1. The annihilation and creation operators will be 
discussed in section |3] in a more general setting. 

3 Other examples of shape invariant "discrete quantum me- 
chanics" 

The other examples of shape invariant "discrete quantum mechanics" are related to a 
two-parameter deformation of the Hermite polynomial, and two- and three-parameter 
deformation of the Laguerre polynomial. The Laguerre polynomials are the eigenfunctions 
of a harmonic oscillator with a centrifugal barrier (1/x 2 potential). They all belong to the 
Askey scheme of hypergeometric orthogonal polynomials 3 . Demonstration of the shape 
invariance and derivation of shift operators and eigenfunctions, etc. go almost parallel with 
the case of the deformed Hermite polynomials in section |2J We discuss them collectively 
by adopting an (almost) self-evident notation. 



Shape Invariant Potentials in "Discrete Quantum Mechanics" 



7 




<t>0 



H a Hi H 2 H 3 

A A + § A+| A+| 

Figure 1. A schematic diagram of the energy levels and the associated Hamiltonian systems to- 
gether with the definition of the A and A* operators and the 'creation' and 'annihilation' operators. 
The parameter set is indicated below each Hamiltonian. 



(i) : two- parameter deformation of the Hermite polynomial, 

a special case of the continuous Hahn polynomial, 

V(x) = V(x; A) = (a + ix)(b + ix), 

4>®{x) = 4>o(x; A) oc |T(a + ix)T(b + ix)\, 

A = (a,6), £ n = £ n (\) = \n(n + 2a + 2b- 1), a,beM 

P n (x) = P( x \x) oc p n (x;a,b,a,b). 

(ii) : two-parameter deformation of the Laguerre polynomial, 
the continuous dual Hahn polynomial, 

(a + ix)(b + ix)(c + ix) 



V(x) = V(x;X) 



(x) = 4> (x; A) oc 



2ix(2ix + 1) 
r(a + ix)T(b + ix)T(c + ix) 



T(2ix) 

\=(a,b, c), £ n = £ n (A) = \n, a,b, c £ M+, 
P n (x) =P^ X) {x) xS n (x 2 ;a,b,c). 
(iii) : three-parameter deformation of the Laguerre polynomial, 
the Wilson polynomial, 

(a + ix)(b + ix)(c + ix)(d + ix) 



V{x) = V{x;\) 



4>o(x) = (f> {x;\) oc 



2ix(2ix + 1) 
r(a + ix)T(b + ix)T(c + ix)T{d 



ix 



T(2ix) 



A = (a, b, c, d), £ n = S n (X) = \n(n + a + b + c + d — 1), a, b,c,d S 
P„(x) = Pi A) (x) oc Ty n (x 2 ; a, 6, c, d). 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



(3.5) 

(3.6) 

(3.7) 
(3.8) 



(3.9) 

(3.10) 

(3.11) 
(3.12) 
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For space reasons we write the ground state wavefunction in terms of the absolute value 
symbol as in Q3.2|) which should read 



|T(a + ix)T(b + ix)\ = \/T(a + ix)T(b + ix)T(a — ix)T(b — ix) 
as in (JTHJ). 

Factorised Hamiltonian Factorisation and consequently the positive semi-definiteness 
of the generic Hamiltonian hold exactly the same as before: 

H = ±^/V(x^e- id *^/W(xj+ ^y/V^xje^y/V^xj- \{V{x) + V*(x)), (3.13) 
H = A^A, (3.14) 

A = A(x; A) d = (e~^ d *y/W(xj- y/V(x)) > ( 3 - 15 ) 

A* =A(x;X)^ = -L^y^)e-i & - y / W(x) et & ). (3.16) 

Here A denotes the (set of) parameters. The ground state 4>q is annihilated by A: 

A(f) O (x) = H<l> (x) = 0, £ = 0. (3.17) 

Verification of the ground state wavefunction (|3.2j) for (i), (|3.6j) for (ii) and (|3.1()j) for (iii) 
is straightforward. The similarity transformed Hamiltonian H = 0q 1 o H o O , (|2.11|) 
determines the other eigenfunctions of the Hamiltonian in the form (j> n (x) oc P n (x)4>o(x), 
HP n = £ n P n ()2.1U|) . They are the difference equation of the form 

V(x)P n (x + V*(x)P n (x + i)- (V(x) + V*{x))P n {x) = 28 n P n {x). (3.18) 

Explicitly they read for the three cases listed above: 

(i) : (a + ix)(b + ix)P n (x — i) + (a — ix)(b — ix)P n {x + i) 

= 2{ab-x 2 +£ n )P n {x). (3.19) 

(ii) : (a + ix)(b + ix)(c + ix)(2ix — \)P n (x — i) + 

(a — ix)(b — ix)(c — ix)(2ix + l)P n (x + i) 

= -2ix [{ab + ac + bc- 2abc + 2(a + b + c)x 2 - x 2 + 2(4x 2 + l)£ n ) P ft (s).(3.20) 

(iii) : (a + ix)(b + ix)(c + ix)(d + ix)(2ix — l)P n (x — i) + 

(a — ix)(b — ix)(c — ix)(d — ix)(2ix + l)P n (x + i) 

= —2ix [(abc + abd + acd + 6ctf — 2abcd + 2(a6 + ac + ad + 6c + bd + cd)x 2 

-(a + b + c + d)x 2 - 2x A + 2(4x 2 + l)£ n ] P n (x). (3.21) 

They admit a degree n polynomial solution in x for (i), and in x 2 for (ii) and (iii). By com- 
paring the coefficients of the leading degree, one obtains the energy eigenvalues £ n (j3.3|) . 
(|3.7jl . (|3.11|) . The solutions form orthogonal polynomials satisfying three term recurrence, 
which will not be shown here, see [3]. 
The corresponding factorisation of H, 

H = BC, C= |(e-s 9 * -e^ 9 *), 5 = -^(a;)^ 9 * _ V*{x)e^ 9x ), (3.22) 
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provides the forward and backward shift operators of the polynomials. The C oper- 
ator, corresponding to the A operator, shifts to the right along the isospectral line: 

C : Pn X \x) — ► P^J^*\x) with | = (5, §,■"■)■ The B operator, corresponding to the 

operator, shifts to the left along the isospectral line: B : P^ l _ l 2 \x) — ► P^\x). For 
each case they are: 

(i) : BCp n (x;\) = £ n p n (x;X) 

<= Cp n (x; A) = (£ n /n)p n -i(x; A + |), Bp n -\{x\ A + |) = np n (x; A), (3.23) 
p n (x + |; A) — p„(x - i; A) = i(n + 2a + 26 - l)p n _i(x; A + §), (3.24) 
(a - ix)(6 - ix)p n _i(x + §; A + \) - (a + ix)(b + ix)p n -i(x - §; A + §) 

= — inp n (x; A). (3.25) 

(ii) : £CS n (x 2 ;A) = ^(x 2 ; A) 

^CS n (x 2 ;A) = -2£ n xS n _i(x 2 ;A+§), -25 x5 n _i(x 2 ; A + §) = S n (x 2 ; A), (3.26) 

S n ((x + §) 2 ; A) - S n ((x - §) 2 ; A) = -2mx5 n _!(x 2 ; A + §), (3.27) 
(a — ix)(6 — ix)(c — ix)S n -i((x + |) 2 ; A + |) 
-(a + ix)(6 + ix)(c + ix)S„_i((x + §) 2 ; A + |) 

= -2ixS n (x 2 ; A). (3.28) 

(iii) : BCW n (x 2 ;X) =£ n W n (x 2 ;X) 

^CW n {x 2 -X) = -28 n xW n ^{x 2 -\+\)i -2BxW n - 1 {x 2 -\+\) = W n {x 2 -\), 

(3.29) 

^ n ((x + §) 2 ;A)-W^((x-§) 2 ;A) 

= -2in(n + a + b + c + d-l) xW n _i(x 2 ; A + J), (3.30) 
(a - ix)(6 — zx)(c - ix)(d — zx)W n _i((x + |) 2 ; A + |) 
-(a + ix)(b + ix)(c + ix)(d + ix)W n - 1 ((x + §) 2 ;A + §) 

= -2ixW„(x 2 ; A). (3.31) 

Shape Invariance For shape invariance, we first need to find the operators Ai, A\ and 
a real constant £ 1 satisfying 

#1 = AA* = Ai + fij (332) 
A 1 = A 1 (x; A) = ±= (e-i d * y/v^xj - e& y/V&j) > (3-33) 
4 = A 1 (x; A)t d ^ f _L(y^) e -^ _ ^i^) e *«°). (3.34) 

In other words, given F(x) = V(a;; A), find a new potential V\(x) = Vi(x; A) satisfying 
V?{x - i)^i(x) = V*(x - §)y(x - §), (3.35) 

Fi(x) + y;(x) = y(x + §) + v*(x - §) + 2^. (3.36) 
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If V\ has the same form as V with a shifted set of parameters A', 

V 1 (x;X) = V(x;X'), (3.37) 

it is shape invariant. Suppose V\ has the form V\(x) = V(x — |)g(x), the above conditions 
(ESP, (EOfl) get slightly simplified: 



#*(x - i)g{x) = 1, 

V(x - i)s(x) + F*(x + f )ff*(x) = F(x + |) + - i) + 2fi. 

The following choice satisfies the above conditions (|3,38|) . (|3.38f) for the three cases (j3.2|) 
(E32J): 

(i) : g(x) = l, £i(X) = a + b, 

V 1 {x;X) = V(x-i;a,b) = V(x;a + \,b + \). 

(h) : 9(x) = 2 -^, £iW = h 

Vi(x; A) = F(x - i;a,b,c)g(x) = V(x;a + -,& + |,c+ 5). 
2ix + 2 

(iii) : g(s) = — , £i(A) = |(a + 6 + c + d), 

A) = V(x - %;a,b,c,d)g(x) = V(x;a + ±, b + \, c + |, d + \). 
With I = (|, |, • • • ), we have collectively 

Vx{x- A) = V(x; A + §), Ax(x; A) = A(x; A + §), tfx(x; A) = H(x; A + §) + £i(A). 

Starting from Vq = V, Hq = H, (f>Q >n = (fr n , let us define V s , H s , ^ Sjn (n > s > 
by step: 

y s+1 (x;A) d =a f y s (x;A + |), 



ff a+ i(x; A) d = A s (x; X)A s (x; A) f + £ S (X), 

<p s+ljn (x; A) d = A s (x; A)0 s , n (x; A), 
f s (A) = f 8 _i(A) + fi(A + (a - l)/2). 

Here A s , A J are defined by 



def 1 

t d -£.f _W / T/ ^. 



A s (x; A) ^Ife-^^/^Aj-e^^A)), 



A s (x; A)t ^ (y^A) e ~5 9 * - ^/(x;A)e3' 

As a consequence of the shape invariance, we obtain for n > s > 0: 
F s (x;A) = l/(x;A+f), 

A s (x; A) = A(x; A + f ), A s (x; A) 1 " = A(x; A + f ) f , 
# s (x; A) = A s (x; X)U s (x; A) + S S (X), 
H s (x; X)(f) s>n {x; A) = £ n (A)0 Sjn (x; A), 
A s (x; A)^ SiS (x; A) = 0, 

A s (x; Xy (j> s +i, n (x; X) = (£„(A) - £ s (A))0 S) „(x; A). 



3.38) 
3.39) 



3.40) 
3.41) 

3.42) 

3.43) 

3.44) 

3.45) 

3.46) 
step 



3.47) 

3.48) 

3.49) 
3.50) 

3.51) 
3.52) 

3.53) 
3.54) 
3.55) 
3.56) 
3.57) 
3.58) 



Shape Invariant Potentials in "Discrete Quantum Mechanics" 



11 



The eigenfunction (f) s ^ n has n — s nodes, or the corresponding polynomial is of degree n — s 
in x for (i), and in x 2 for (ii) and (hi). In the latter case, the zeros on a half line, say 
x > 0, count as the nodes. 

From these we obtain formulas 

(f> s>n (x;X) = As^x; X) ■ ■ ■ A 1 (x; \)A (x; \)(j) n (x; X), (3.59) 

4)(s;A)t Ai(a;A)t A^faA^ 
^ A) = £ n {X)-£, { X) £ n (X) - £\{X) ■ ■ ■ £ n (X)-£ n ^(X) ^ ^ 

corresponding to (|2.34j) and (|2.42j) discussed in section [2j The former (j3,59|) gives the 
eigenfunction cj) s ^ n of the s-th Hamiltonian H s along the isospectral line with energy £ n , 
starting from <p n of the original Hamiltonian H by repeated application of the A oper- 
ators. The latter (|3.6()fl . on the other hand, expresses the n-th eigenfunction <p n of the 
original Hamiltonian, starting from the explicitly known ground state 4> njn of the n-th 
Hamiltonian H n by repeated application of the A^ operators. The latter formula (|3.6()[) . a 
simple generalisation of the well-known formula for the harmonic oscillator |n) tx (a^)™|0), 
could also be understood as the generic form of the Rodrigue's formula for the orthogonal 
polynomials. 

In order to define the annihilation and creation operators, let us introduce normalised 
basis {<t> s ,n}n>s for each Hamiltonian H s . Ordinarily, the phase of each element of an 
orthonormal basis could be completely arbitrary. In the present case, however, the eigen- 
functions are orthogonal polynomials. That is, they are real and the relations among 
different degree members are governed by the three term recurrence relations. So the 
phases of {<p s ,n}n>s are fixed. Let us introduce a unitary (in fact an orthogonal) operator 
U s mapping the s-th orthonormal basis {4> s , n }n>s to the (s + l)-th {(p s+ i tn } n > s+ i (see Fig. 
1 and for example [TBI IT7] ): 

U s 4> Stn = s+ l >n+ l, U^ s+ l in+ l = (j) s , n . (3.61) 

We denote that Uq = U. Roughly speaking U increases the parameters from A to A + ~: 



Let us introduce an annihilation a and a creation operator for the Hamiltonian H as 
follows: 

a = a(x; A) = f ll' f A(x; A), a) = a(x; A)* = f A(x; X^U. (3.62) 
It is straightforward to derive 

H = a f a, (3.63) 
[a,a^]4> n (x; A) = (£„+i(A) - £ n {X))4> n {x; A). (3.64) 

For the harmonic oscillator (|2.1|) U =id. and we recover the known result. For the linear 
spectrum £ n = n,n/2 of the deformed Hermite polynomial (|2.3f) and the two parame- 
ter deformation of the Laguerre polynomial (|3.5|) - (|3.8|) . a and a) have the same (up to 
rescaling) commutation relations as those of the harmonic oscillator. Essentially the same 
arguments and results hold for the annihilation a s = u\A s (x\ A) and the creation operator 
at = AAx\ AW, for the Hamiltonian H q . 
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4 Dynamical Background 

Here we will show briefly a logical (dynamical) path that led to the shape invariant differ- 
ence equations introduced in section [21 The equilibrium position of an n-particle A-type 
Calogero system is determined by 

n 1 

=Xj, j = l,...,n, (4.1) 

nr. a — nri. J 



fc=i J 
k±j 

after adjustment of the coupling constants and rescaling of the variables [SJ E|. They 

describe the zeros of the Hermite polynomial, since H n (x) = 2 n Y\™ = ±(x — Xj) satisfies the 
differential equation 

H„(x) - 2xH' n {x) + 2nH n (x) = (4.2) 

and the three term recurrence j^j n > 0, 

H n+1 (x)-2xH n (x) + 2nH n ^(x) = 0, H = 1, = 0. (4.3) 

The equilibrium position of an n-particle A-type van Diejen system (or the rational R-S 
system with a linear confining potential) is determined by 

nXj-x k -iV5 l-iy/8xj . 
— ~r = - ~Jt~ > J = !,■■■,«. (4-4) 
fc=1 £j — Xk + w o 1 + ivd £j 

after adjustment of the coupling constants and rescaling of the variables |llj . The corre- 

Ii 



sponding deformed polynomial H n (x, 5) d = 2 n YYj=i( x ~ x j) satisfies the difference equation 



2i 

and the three term recurrence n > 0, 



x+ -^)H n {x + iV8,8) - (x - -^jH n (x -iy/5,5) = -^(1 + n6)H n (x, S), (4.5) 



H n+1 (x,5)-2xH n {x,5) + (2n + n{n-l)5)H n ^ l (x,5) = 0, H = 1, = 0. (4.6) 

It is elementary to show that (|4~ij) . ijOjl . reduce to (|4~T|) . (j4~2~|) . in the zero de- 
formation limit 5 — > 0; lim^o H n (x, 5) = H n (x). The relationship between the deformed 
Hermite polynomial and the special case of the Meixner-Pollaczek polynomial discussed 
in section |2] is 

H n (x,5)=n\V6 n P { J\^=;^). (4.7) 

It is instructive to write down the classical Hamiltonian of the n-particle A-type van Diejen 
system 

H(p,q) = £ [ cosh Pi y/Vj(q)Vf(q) ~ ^(Vj(q) + V*(q))j , (4.8) 



o I t=i \ H ~ Ik 
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in which a and g are coupling constants. In fact, the Hamiltonian (|2.3j) is the single 
particle (n = 1) case with pj = —id/dqj and proper ordering. We refer to for similar 
orientation of the dynamical systems introduced in sectional 



5 Comments and Discussion 

Several examples of shape invariant difference equations are discussed in some detail. They 
are related to deformation of the harmonic oscillator without /with a centrifugal potential 
and their eigenfunctions are deformed Hermite and Laguerre polynomials belonging to 
the Askey-scheme of hypergeometric orthogonal polynomials. They arise in the problems 
of determining the equilibrium positions of the Ruijsenaars-Schneider-van Diejen systems, 
which are integrable deformation of the celebrated Calogero-Sutherland systems of exactly 
solvable multi-particle quantum mechanics. Here we have treated rational potentials V(x) 
only. Obviously the method works for a wider range of potentials, the trigonometric, 
hyperbolic, elliptic, etc. We have not discussed the shape invariance of the discrete systems 
corresponding to the equilibria of the trigonometric Ruijsenaars- Schneider systems. These 
have various kinds of deformed Jacobi polynomials as eigenfunctions We will report 
on this subject elsewhere. 

In the literature, shape invariance is discussed almost always within the context of 
'supersymmetric quantum mechanics '[.1.. We presume that this might be the psychological 
barrier for considering the shape invariance in discrete quantum mechanics, or shape 
invariant difference equations. As is well-known, supersymmetry is the 'square root' of the 
space-time symmetry, say the Poincare symmetry, which is usually lost if the space-time 

is discretised In fact, it is important to realise that main results of 'supersymmetric 

quantum mechanics' are already contained in Crum's theorem without supersymmetry or 
shape invariance. Among them are: the existence of associated isospectral Hamiltonians 
H\, . . . , H s , . . . , as many as the number of discrete levels of Hq and the formulas of their 
eigenfunctions. For example, <p s ^ n of H s with the eigenvalue S n is expressed in terms of 
the eigenfunctions tfio, <j>i,. . . , <p s -i, <j) n of Hq (Fig. 1): 

4> s ,n = -r^f 1 , n > s > 1, (5.1) 
W s = W[(t> ,...,(f> s -i], W.,„ = W[^o,...,^_ 1 ,0 n ]. (5.2) 

Here the Wronskian W is defined as usual W[/i,...,/ n ] = det ( /■* *M . This 

V J J l<i,j<n 

formula, translated in the discrete dynamics context, corresponds to (|3.59[) . 

We do believe dynamical (Hamiltonian) interpretation of discrete quantum systems 
(difference equations) would be useful and fruitful. One of our starting points, equations 
determining a certain equilibrium, eg. (|4.4jl . are called Bethe ansatz like equations. The 
present problems could be considered as those of finding associated polynomial solutions 
of Bethe ansatz like equations. They occur in many branches of theoretical physics, for 
example, the quasi-exactly solvable single and multi-particle quantum systems ^H] on top 
of the well-known integrable spin chains |19j . 
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